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HYPERSONIC  AERODYNAMICS  FELLOWSHIPS 
FINAL  REPORT 

by 

John  D.  Anderson,  Jr.,  Professor 
Department  of  Aerospace  Engineering 
University  of  Maryland 
College  Park,  Maryland  20742 


m 


This  is  a  final  report  on  ARO  Grant  DAAL03-88-G-0083.  This  research 
program  was  to  provide  graduate  fellowships  to  students  in  the  area  of 
hypersonic  aerodynamics.  Indeed,  the  Army  Research  Office  can  take  some 
pride  in  the  fact  that  this  funding  was  the  leading  edge  of  a  renewed  emphasis 
in  hypersonic  aerodynamics  in  universities  during  the  middle  1980's,  at  a 
time  when  this  discipline  was  seriously  languishing.  In  particular,  this  grant 
served  to  spur  the  growth  of  the  Hypersonics  Group  at  the  University  of 
Maryland-  a  group  which  now  numbers  17  full-time  graduate  students.  In 
this  sense,  we  wish  to  express  our  gratitude  to  ARO  for  helping  our  graduate 
program  in  hypersonic  aerodynamics  in  such  a  pivotal  manner. 

The  following  students  were  funded  under  the  ARO  Fellowship 
Program  at  Maryland.  Also  listed  are  their  current  status. 


1.  Dr.  Kevin  Bowcutt.  Received  his  Ph.D.  under  the  ARO  Fellowship 
Program,  and  is  now  a  lead  research  engineer  with  Rockwell  in 
their  National  Aerospace  Plane  Program.  Dr.  Bowcutt  was  selected 
as  the  Rockwell  Engineer  of  the  Year  in  1990  —  a  very  prestigious 
award  for  such  a  young  engineer. 
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2.  Dr.  Griffin  Corpening.  Received  his  Ph.D.  under  the  ARO 
Fellowship  Program,  and  is  now  a  research  leader  at  the  Johns 
Hopkins  Applied  Physics  Laboratory  (APL)  in  their  hypersonic 
SC  RAM  JET  program.  Last  year,  Dr.  Corpening  was  awarded  APL's 
Best  Paper  of  the  Year  Award,  for  an  AIAA  paper  based  on  his  Ph.D. 
dissertation.  (This  paper  is  included  as  Appendix  B  in  this  Final 
Report.) 

3.  Mr.  Scott  Seifert.  Finished  all  requirements  for  his  Ph.D.  except  for 
his  Oral  Comprehensive  Examination  and  his  Dissertation.  Mr. 
Seifert  is  now  continuing  to  work  towards  finishing  these 
requirements,  while  at  the  same  time  accepting  employment  at  the 
Applied  Physics  Laboratory. 

The  following  accomplishments  have  been  made  as  a  direct  result  of 
this  ARO  Fellowship  Grant: 

1.  Two  Ph.D.'s  and  one  near  Ph.D  have  been  added  to  the  United 
State's  engineering  work  force  in  hypersonics,  at  a  critical  time 
when  they  were  badly  needed.  Moreover,  these  graduates  are  top- 
notch  people  —  the  creme  of  the  crop.  Both  Drs.  Bowcutt  and 
Corpening  graduated  with  perfect  4.0  Grade  Point  Averages  — 
straight  A  s. 

2.  Dr.  Bowcutt's  Ph.D.  research  dealt  with  a  new  approach  to 
hypersonic  waverider  design,  and  led  to  a  new,  unique  family  of 
hypersonic  waveriders  called  viscous  optimized  waveriders.  These 
University  of  Maryland  waveriders  have  prompted  a  new  surge  in 
waverider  activity  around  the  country.  Our  waverider  computer 
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program  (MAXWARP)  is  now  being  used  by  15  companies  and 
government  laboratories.  Two  University  of  Maryland  waveriders 
have  been  successfully  tested  in  the  NASA  Langley  hypersonic 
wind  tunnels,  validating  our  analysis  and  multiplying  interest  in 
these  waveriders.  Langley  is  now  embarking  on  their  own  high- 
priority  waverider  program,  as  a  result.  In  October  1990,  the 
University  of  Maryland  hosted  the  1st  International  Hypersonic 
Waverider  Symposium,  which  was  an  acclaimed  success.  There 
have  been  three  articles  in  Aviation  Week  over  the  past  four 
months  describing  the  University  of  Maryland  Waverider  Program. 
All  of  this  was  started  with  the  ARO  Fellowship  Grant.  We 
sincerely  believe  that  ARO  can  take  pride  in  this  accomplishment  — 
the  ARO  Grant  was  a  seed  upon  which  a  major  program  has  been 
built  -  a  program  with  far-reaching  results  affecting  many  people  in 
our  national  hypersonics  effort. 

3.  Dr.  Corpening's  Ph.D.  research  was  an  application  of  computational 
fluid  dynamics  to  a  flow  problem  of  interest.  The  flow  problem  was 
the  interaction  of  a  vortex  with  an  oblique  shock  wave  at 
hypersonic  speeds.  Dr.Corpening  wrote  his  own  three-dimensional, 
finite-volume  computer  program  to  solve  this  problem.  The  results 
indicated  a  major  interaction  between  the  vortex  and  the  shock,  the 
first  such  data  at  hypersonic  speeds. 

The  details  of  this  technical  work  are  contained  in  Appendices  A  and  B. 
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These  are: 


Appendix  A:  Bowcutt,  K.G.,  Anderson,  J.D.  and  Capriotti,  D.,  "Viscous 
Optimized  Hypersonic  Waveriders,"  AIAA  Paper  No.  87-0272. 

Appendix  B:  Corpening,  G.  and  AndersonJ.,  "Numerical  Solutions  to 
Three-Dimensional  Shock  Wave/ Vortex  Interaction  at  Hypersonic 
Speeds,"  AIAA  Paper  No.  89-0674. 

Hence,  no  further  elaboration  will  be  given  here.  The  remainder  of  this 
report  is  represented  by  the  technical  descriptions  given  in  Appendices  A  and 
B. 
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APPENDIX  A 


AIAA-87-0272 

Viscous  Optimized  Hypersonic 
Waveriders 

K.  G.  Bowcutt,  J.  D.  Anderson  and 
D.  Capriotti,  Univ.  of  Maryland, 
College  Park,  MD 


AIAA  25th  Aerospace  Sciences  Meeting 

January  12*15, 1987/Reno,  Nevada 


For  permission  to  copy  or  republic*,  contact  the  American  InsMuto  of  Aoroaautict  and  Aatronautlca 

1193  Broadway,  Now  York,  NY  10019 


VISCOUS  OPTIMIZE!!  HYPERSONIC  MAVERIDERS 


uy 


Kevin  G.  Bowcutt,*  John  0.  Anderson,  Jr.,**  and  Oiego  Caprlotti*** 


Department  of  Aerospace  Engineering 
University  of  Maryland 
College  Park,  Maryland  20742 


"I  t  is  neatly  not  iorcseeable  that  an  '  optimized’  ca Iculattd  a  hope,  could  do  anything 
mo> te  than  give  a  guide  to  the  designer.  However,  it  is  only  a  guide  that  the  designer 
actually  i amts  &rom  the  aerodynamicist.  He  would  xeally  be  a  little  embarrassed  to  be 
o^ered  a  per&eci  aerodynamic  shape,  which  he  would  then  have  to  carve  holer  in,  add 
gainings,  and  so  on,  in  older  to  satisfy  iuch  mundane  requirements  as  that  the  pilot 
should  be  able  to  see  where  he  is  going  ox  that  people  have  somewhere  convenient  to  get 


in  and  out.” 


ABSTRACT 


A  family  of  optimized  hypersonic  wavertders 
Is  generated  and  studied  wherein  detailed  viscous 
effects  are  Included  within  the  optimization  pro¬ 
cess  Itself,  This  Is  In  contrast  to  previous 
optimized  waverlder  work,  wherein  purely  Invlscld 
flow  Is  used  to  obtain  the  waverlder  shapes.  For 
the  present  waverlders,  the  undersurface  Is  a 
streamsurface  of  an  Invlscld  conical  flowfleld, 
the  upper  surface  Is  a  streamsurface  of  the 
Invlscld  flow  over  a  tapered  cylinder  (calculated 
by  the  axl symmetric  method  of  characteristics), 
and  the  viscous  effects  are  treated  by  Integral 
solutions  of  the  boundary  layer  equations. 
Transition  from  laminar  to  turbulent  flow  Is 
Included  within  the  viscous  calculations.  The 
optimization  Is  carried  out  using  a  non-linear 
simplex  method.  The  resulting  family  of  viscous 
hypersonic  waverlders  yields  predicted  high 
values  of  lift/drag,  high  enough  to  break  the 
"1/0  barrier"  based  on  experience  with  other 
hypersonic  configurations.  Moreover,  the 
numerical  optimization  process  for  tlu  viscous 
waverlders  results  In  distinctly  different 
shapes  compared  to  previous  work  with  Invlscld 
-designed  waverlders.  Also,  the  fine  details 
of  the  viscous  solution,  such  as  how  the 

*  Graduate  hypersonic  Aerodynamics  Fellow, 
sponsored  by  the  Array  Research  Office, 

Or.  Robert  Singleton,  monitor.  Member,  AIAA. 
Presently,  Engineer,  Rockwell  International, 
Los  Angeles. 

**  Professor.  Fellow,  AIAA.  Currently  on 
sabbatical  as  the  Charles  Lindbergh 
Professor  at  the  National  Air  and  Space 
Museum,  Smithsonian  Institution. 

***  Graduate  Research  Assistant.  Student  member, 
AIAA. 


P.L.  Rpe 

Royal  Aircraft  Establishment 
January  1970 


shear  stress  Is  distributed  over  the  surface,  and 
the  location  of  transition,  are  crucial  to  the 
details  of  the  resulting  waverlder  geometry. 

I.  INTRODUCTION 

Over  the  past  few  years.  Interest  In  all 
aspects  of  hypersonic  flight  has  grown  explosi¬ 
vely,  driven  by  new  vehicle  concepts  such  as  the 
National  Aerospace  Plane  (HASP),  aero-asslsted 
orbital  transfer  vehicles  (AOTV's),  the  hyper¬ 
sonic  transport  (the  "Orient  Express"),  and 
hypersonic  missiles,  to  name  just  a  few.  An 
extended  discussion  of  these  concepts,  as  well  as 
a  survey  of  hypersonic  aerodynamic  research 
contrasting  the  "old"  with  the  "new*  hypersonlcs, 
is  given  In  Ref.l.  Hence  no  further  elaboration 
will  be  given  here. 

The  present  paper  deals  with  a  class  of 
advanced  hypersonic  lifting  configurations.  To 
help  understand  the  motivation  for  the  present 
work,  the  following  background  Is  given.  For  a 
lifting  aerodynamic  body.  It  Is  well-known  that 
high  maximum  llft-to-drag  ratios,  (L/D)max«  ar# 
very  difficult  to  obtain  at  hypersonic  speeds, 
due  to  the  presence  of  strong  shock  waves  (hence 
high  wave  drag)  and  massive  viscous  effects.  At 
supersonic  and  hypersonic  speeds,  the  most  effj- 
clent  lifting  surface  Is  the  Infinitely  thin  flat 


pi «t*;  the  Invlscld  hypersonic  torodynMlc  pro¬ 
perties  of  a  flat  plate  art  shown  as  th*  solid 
curves  In  Flg.l,  based  on  the  Newtonian  Unit  of 
free  stream  Mach  number  and  y  ■  Cn/Cy»l. 

Note  that  I/O  theoretically  approaches  infinity 
as  the  a.-.gle-of -attack,  a,  approaches  zero.  In 
reality,  viscous  effects  will  cause  I/O  to  peak 
at  low  values  of  <1,  and  to  go  to  zero  as  a-*0. 

This  Is  Illustrated  by  the  dashed  line  In  Flg.l, 
which  shows  the  variation  of  I/O  modified  by  skin 
friction  as  predicted  by  a  reference  temperature 
method  described  in  Ref.  2.  Although  the  Inf  Ini t- 
1y  thin  flat  plate  shown  In  Flg.l  Is  the  most 
effective  lifting  surface  aerodynamical ly.  It  ts 
the  least  effective  In  terms  of  volume  capactty. 

It  goes  without  saying  that  all  practical  flight 
vehicles  must  have  a  finite  volume  to  carry  fuel, 
payload,  etc.  Hence,  th*  flat  plate  results,  al¬ 
though  instructive,  are  of  academic  Interest 
only.  In  contrast.  Fig.  2  shows  values  of 
(l/D)max  v*r*“*  the  volume  parameter  Vz/1 2 3/S  for 
several  generic  hypersonic  configurations,  obtaln- 
from  Ref.  3.  Here,  V  Is  the  body  volume  and  S  is 
the  planfone  area.  Note  from  Fig. 2  that  typical 
hypersonic  values  of  (l/0)MX  rang*  from  4  to  6 
for  such  lifting  bodies  at  th*  conditions  shown. 
These  values  are  also  typical  of  the  hypersonic 
transport  conf Iguratlon  studied  In  Ref.  4.  Clear¬ 
ly,  values  of  (L/D)Max  for  hypersonic  vehicles 
are  substantially  lower  than  those  for  conven- 
al  subsonic  and  low  supersonic  airplanes.  (For 
example,  the  maximum  I/O  values  for  the  World  War 
II  8o*1 ng  B-29  and  th*  contemporary  General  Dy¬ 
namics  F-IU  are  lb.8  and  15.8  respectively,  as 
obtained  from  Ref.  5).  Indeed,  as  N.  Increases 
across  the  supersonic  and  hypersonic  regimes, 
there  Is  a  general  empirical  correlation  for 
(L/0)fflax  based  on  actual  flight  vehicle  experi¬ 
ence,  given  by  Kuchemann*: 


4(H-*3) 

C-^max  "  IT 

This  variation  Is  shown  as  the  solid  curve  In 
Fig. 3.  Also  shown  are  a  number  of  data  points 
for  various  previous  hypersonic  vehicle  configu¬ 
rations  at  various  Reynolds  numbers  (the  open  sym¬ 
bols),  as  well  as  new  results  from  the  present  in¬ 
vestigation  (the  solid  symbols).  Fig.  3  is  pivo¬ 
tal  to  the  present  paper,  and  will  be 
discussed  at  length  In  subsequent  sections.  How¬ 
ever,  at  this  stage  In  our  discussion.  Fig. 3  is 
used  to  Illustrate  only  th*  following  aspects: 

1.  The  solid  curve  represents  a  type  of  "l/D 
barrier"  for  conventional  hypersonic  vehi¬ 
cles,  which  Is  difficult  to  break. 

2.  OaM  for  conventional  hypersonic  vehicles, 
shown  as  th*  open  circles, 

form  an  almost  random  "shotgun"  pattern 
which,  for  the  most  part,  falls  below  th* 
solid  curve. 

(The  numbers  adjacent  to  these  open 
circles  pertain  to  specific  reference  num¬ 
bers  Itemized  In  Ref.  7,  which  should  be 
consulted  for  details.) 

3.  The  solid  symbols  pertain  to  th*  present 
study,  and  represent  a  new  class  of  hyper¬ 
sonic  configurations  which  break  th*  "l/D 
barrier."  These  configurations  are  conical 
flow  waverlders  that  are  optimized  with  de¬ 
tailed  viscous  effects  Included  directly  In 
the  optimization  process. 


To  help  understand  the  contribution  mad*  by 
th*  present  work,  let  us  briefly  review  the  gener¬ 
al  concept  of  waverlders.  if.  1959.  the  design  of 
three-dimensional  hypersonic  vehicles  which  support 
planar  attajh|d  shock  waves  was  Introduced  by 
Nonwel  Her,  ’  who  hypothesized  that  streamsurfaces 
from  the  flow  behind  a  planar  oblique  shork  could 
be  used  as  supersonic  lifting  surfaces.  This  led 
to  a  class  of  vehicles  with  a  caret-shaped  trans¬ 
verse  cross-section  and  a  delta  planform—  the  so- 
called  caret  wing  as  shown  In  Fig. 4.  Here,  the 
body  surface  Is  generated  by  stream  surfaces  behind 
a  planar  oblique  shock  wave.  The  shock  wave  Is 
attached  to  the  sharp  leading  edges  at  the  design 
Hach  number,  and  hence  no  flow  spillage  takes 
place  around  the  leading  edge.  The  lift  is  high 
due  to  the  high  pressures  behind  a  two-dimensional 
planar  shock  wave,  exerted  on  the  lower  surface  of 
the  vehicle.  Because  the  body  appears  to  be  rid¬ 
ing  on  top  of  the  attached  shock  wave.  It  Is 
called  a  "waverlder".  The  aerodynamic  advantages 
of  such  waverlders  are  listed  In  Ref.  1,  and  are 
discussed  in  great  detail  In  Refs.  7  and  10.  In 
short,  without  repeating  the  details  here,  at  a 
given  lift  coefficfent,  caret  waverlders  theore¬ 
tically  operate  at  higher  l/D  values  than  other 
hypersonic  configurations. 

Expanding  on  this  philosophy,  other  types  of 
flowfields  can  be  used  to  generate  waverlders.  For 
example,  any  streamsurface  from  the  supersonic 
flow  over  an  axl symmetric  body  can  be  used  to 
generate  a  waverlder  with  an  attached  shock  wav* 
along  its  complete  leading  edge.  Work  on  such 
waverlders  was  first  carried  out  in  Britain,  as 
nicely  summarized  In  Ref.  11,  where  th*  flow  over 
a  right-circular  con*  at  zero  degrees  angle  of 
attack  Is  used  to  generate  a  class  of  "conical 
flow*  waverlders.  Still  later,  waverlders  were 
generated  from  Inclined  cfrcular  and  elliptic 
cones,  and  axlsyimnetrlc  bodies  with  longitudinal 
curvature  by  Rasmussen  et  al.  *  ,  using  hyper¬ 

sonic  small  disturbance  theory.  This  work  was 
further  embellShed  by  the  search  for  optimize 
waverlder  shapes.  For  example  Cole  and  Zien 
found  optimized  waverlders  derived  from  axlsym- 
metrlc  bodies  with  longitudinal  curvature  by  using 
hypersonic  small  disturbance  theory  to  generate 
Invlscld  flow  solutions,  and  then  utilizing  the 
calculus  of  variations  to  obtain  the  optimum 
waverlder  shapes.  Later,  Kim  et  al.  used  the  same 
philosophy  to  derive  optimum  waver1derst Jrom 
flowfields  about  unyawed  circular  cones  ,  and 
yawed  circular  and  elliptic  cones  .  In  Refs.  14- 
16,  the  advantage  of  hypersonic  small  disturbance 
theory  Is  that  analytic  expressions  are  obtained 
for  surface  pressure  distributions,  hence  lift  and 
wave  drag,  thus  enabling  the  application  of  th* 
calculus  of  variations  for  optimization. 

Unfortunately,  to  date  the  potential  superiority 
of  waverlders  as  hypersonic  high  L/D  shapes  has 
not  been  fully  demonstrated,  either  in  the  wind 
tunnel  or  in  flight.  A  basic  problem  arises 
because  of  the  tendency  for  waverlders  to  have 
large  wetted  surface  areas,  which  leads  to  large 
friction  drag.  411  previous  waverlder  optimiza¬ 
tion  work  (such  as  Refs.  14-16)  has  been  based  on 
the  assumption  of  invlscld  flow,  after  which  an 
estimate  of  skin  friction  for  the  resulting  con¬ 
figuration  Is  sometimes  added.  As  a  result,  the 
real  aerodynamic  performance  of  th*  resulting  opti¬ 
mum  configuration  usually  falls  short  of  Its  expec¬ 
tations. 
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The  purpose  of  the  present  work  Is  to  remove 
this  deficiency.  In  particular,  a  series  of  con¬ 
ical-flow  generated  waverlders  are  optimized  for 
maximum  I/O  wherein  detailed  viscous  effects  (In¬ 
cluding  boundary  layer  transition)  are  Included 
within  tne  optimization  process  itself.  This 
leads  to  a  new  class  of  waverlders  where  the  opti¬ 
mization  process  Is  trying  to  reduce  the  wetted 
surface  area,  hence  reducing  skin  friction  drag, 
while  at  the  same  time  maximizing  l/D.  Because 
detailed  viscous  effects  can  not  be  couched  in 
simple  analytical  forms,  the  formal  optimization 
methods  based  on  the  calculus  of  variations  can 
not  be  used.  Instead,  in  the  present  work  a 
numerical  optimization  technique  is  u«d,  based  on 
the  simplex  method  by  Nelder  and  Head1  .  By  using 
a  numerical  optimization  technique,  other  real 
configuration  aspects  can  be  included  in  the  ana* 
lysis  in  addition  to  viscous  effects,  such  as 
blunted  leading  edges,  and  an  expansion  upper  sur¬ 
face  (In  contrast  to  the  standard  assumption  of  a 
free  stream  upper  surface,  l.e.,  an  upper  surface 
with  all  generators  parallel  to  the  freestream 
direction).  The  results  of  the  present  study  lead 
to  a  new  class  of  waverlders,  namely  "viscous 
optimized"  waverlders.  Moreover,  these  waverlders 
appear  to  produce  relatively  high  values  of  (L/0), 
as  will  be  discussed  in  subsequent  sections. 

II.  analysis 


For  the  present  waverlder  configurations,  the 
following  philosophy  Is  followed: 

1.  The  lower  (compression)  surface  is  generated 
by  a  streamsurface  behind  a  conical  shock 
wave.  The  Invlscld  conical  flowfleld  Is 
obtained  from  the  numerical  solution  of  the 
Taylor-Maccoll  equation,  derived  for  example 
In  Ref.  18. 

2.  The  upper  surface  Is  treated  as  an  expansion 
surface,  generated  In  a  similar  manner  from 
the  Invlscld  flow  about  a  tapered,  axisym- 
metric  cylinder  at  zero  angle  of  attack,  and 
calculated  by  means  of  the  axlsymmetrlc 
method  of  characteristics. 

3.  The  viscous  effects  are  calculated  by  means 
of  an  integral  boundary  layer  analysis 
following  surface  streamlines.  Including 
transition  from  laminar  to  turbulent  flow. 

A.  Blunt  leading  edges  are  Included  to  the 
extent  of  determining  the  maximum  leading 
edge  radius  required  to  yield  acceptable 
leading  edge  surface  temperatures,  and  then 
estimating  the  leading  edge  drag  by  modified 
Newtonian  theory. 

5.  The  final  waverlder  configuration,  optimized 
for  maximum  L/0  at  a  given  Mach  number  and 
Reynolds  number  with  body  fineness  ratio  as 
a  constraint.  Is  obtained  from  the  numerical 
simplex  method  taking  into  account  all  the 
effects  Itemized  In  1-4  above  within  the 
optimization  process  itself. 

The  following  paragraphs  describe  each  of  the 
above  Items  In  more  depth;  for  a  highly  detailed 
discussion,  see  Ref.  7. 


A.  Invlscld  Flow  --  Lower  Surface 


The  waverlder's  lower  surface  is  generated 
from  a  streamsurface  behind  a  conical  shock  wave 
supported  by  a  hypothetical  right  circular  cone 
at  zero  angle  of  attack.  The  hypothetical  cone 
and  its  flowfleld  Is  shown  In  Fig.  s,  where  ec  is 
tne  cone  semi-angle  and  os  is  the  wave  angle.  The 
inviscid  conical  flow  is  obtained  from  the  Taylor. 
Haccoll  equation1® 
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solved  by  a  standard  Runge-Kutta,  forth-order 
accurate  numerical  technique,  namely  the  RKF45 
algorithm  obtained  from  Ref.  19.  In  Eq.(l),  Vp  Is 
the  component  of  nondimensional  flow  velocity 
•long  a  conical  ray,  a  is  the  angle  of  the  ray 
referred  to  the  cone  axis,  and  y  1$  the  ratio  of 
specific  heats. 


Any  streamsurface  from  this  flowfleld  can 
represent  the  wing  undersurface  of  a  waverlder, 
as  shown  In  Fig.  6.  (For  purposes  of  Illustra¬ 
tion,  Fig.  fi  also  shows  the  waverlder  upper  sur¬ 
face  as  a  freestream  surface,  but  this  Is  only 
one  of  many  possible  choices  for  the  upper 
surface.)  Any  particular  undersurface  Is  uni¬ 
quely  defined  by  the  intersection  of  the  conical 
flow  streamsurface  with  the  conical  shock  wave, 
as  shown  by  the  curve  labeled  "leading  edge"  In 
Fig.  7.  Let  us  examine  Mg. 7  more  closely.  It 
is  a  front- view  of  the  hypothetical  conical 
flowfleld.  Illustrating  the  cone  apex  at  the 
center,  and  both  the  cone  base  and  shock  base  at 
some  arbitrary  distance  downstream  of  the  apex. 
Consider  a  curve  In  this  front-view,  lying  below 
the  apex  (or  even  Including  the  apex),  as  shown 
by  the  curve  labeled  "leading  edge*.  Row 
construe  this  curve  as  a  trace  on  the  conical 
shock  wave  Itself,  and  visualize  streamlines 
trailing  downstream  from  this  trace;  the 
resulting  streamsurface  Is  the  waverlder  under, 
surface  sketched  in  Fig.  6.  Indeed,  the  curve 
labeled  "leading  edge"  In  Fig.  7  Is  simply  the 
forward  projection  of  the  waverlder  leading  edge 
on  the  cross-flow  (x-y)  plane.  This  curve  Is 
treated  as  completely  general,  except  for  the 
constraints  that  It  be  synmetric  about  the  y-z 
plane,  and  that  It  lie  entirely  below  the  x-z 
plane  to  ensure  that  the  waverlder  undersurface 
Is  a  compression  surface.  Also  In  Fig.  7,  note 
the  curve  labeled  "trailing  edge*.  This  Is  the 
Intersection  of  the  particular  conical  flow 
streamsurface  with  the  plane  of  the  shock  base, 
and  It  represents  the  bottom  surface  trailing 
edge  of  the  waverlder.  This  Is  the  shape  of  the 
bottom  of  the  waverlder  base,  as  sketched  In  Fig. 
rr  Returning  to  Fig. 7,  the  area  between  the 
"leading  edge*  and  "trailing  edge*  curves  Is  the 
forward  projection  of  the  entire  waverlder 
compression  surface.  Moreover,  the  dashed  line 
emanating  from  the  cone  apex  In  Fig.  7  Is  the 
forward  projection  of  a  conical  flow  streamline; 
hence,  that  portion  of  the  dashed  line  contained 
between  the  "leading  edge"  and  "trailing  edge" 
curves  is  the  projection  of  a  particular 
streamline  along  the  waverlder  undersurface,  from 
the  leading  edge  to  the  trailing  edge. 
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8.  Inviseld  flow  -«  Upper  Surface 

In  w*t  previous  waverider  work,  the  upper 
surf ic“  's  treated  as  a  freestream  surface,  as 
tllustratert  n  Fig.  6.  Here,  the  upper  surface 
pressure  Is  freestream  pressure,  p..  However,  If 
tne  upper  surface  Is  made  an  expansion  surface, 
where  p<p..  then  a  small  but  meaningful  contribu¬ 
tion  to  L/O  can  be  obtained.  This  approach  is 
taken  here.  Similar  to  the  philosophy  for  the 
lower  surface,  the  upper  surface  Is  a  streamsur- 
face  “carved"  from  a  known  expansion  flow.  The 
hypothec i cal  expansion  body  chosen  here  is  a  cir¬ 
cular  cylinder  of  given  radius;  the  cycllnder  is 
aligned  parallel  to  tne  flow  and,  at  some  point, 
is  tapered  parabolically  to  a  smaller  radius.  The 
result  is  an  ax  i  symmetric  expansion  flow,  where 
the  domain  of  expansion  is  bounded  by  a  freestream 
Mach  cone  centered  on  the  cy  Under  axis,  as  shown 
in  Fig.  8.  Parabolic  taper  was  chosen  because  it 
is  relatively  simple,  and  the  resulting  expansion 
body  slope  is  everywhere  continuous.  Once  the 
expansion  body  is  chosen,  it  remains  only  to 
geometrical ly  position  the  expansion  region  rela¬ 
tive  to  the  lower  surface,  choose  the  initial  and 
final  cylinder  radii,  solve  the  Inviseld  expansion 
flow,  then  cut  a  streamsurface  from  that  f  m  to 
serve  as  the  waverider  upper  surface.  This  basic 
idea  was  first  developed  for  two-dimensional 
expansions  by  Flower*  ,  and  later  for  uisymmetric 
expansions  by  Moore  . 

The  axi symmetric  flow  is  calculated  from  the 
axisymnetric  method  of  characteristics,  using  the 
two-step  predictor-corrector  iteration  scheme  of 
Ferri  .  The  details  involving  the  matching  of 
the  resulting  expansion  surface  with  the  conical 
flow  compression  surface  are  straight  forward,  hut 
lengthy.  Considering  that  the  expansion  surface 
contrioutes  only  about  10%  to  the  value  of  (I/O), 
no  further  space  for  Its  discussion  is  justified 
In  the  present  paper;  for  the  complete 
discussion,  see  Ref.  7. 

C.  Leading-Edge  Bluntness 

Waveriders,  by  design,  have  sharp  leading 
edges  that  support  attached  shock  waves. 

However,  for  flight  Mach  numbers  above  five,  the 
temperatures  for  sharp  leading  edges  will  exceed 

practical  limits  of  most  structural 
materials.  Th1>  leads  to  the  need  for  blunt 
leading  edges  with  sufficiently  large  radii  such 
that  the  aerodynamic  heat  flux  is  reduced  to 
reasonable  levels.  However,  at  the  sane  time  the 
leading  edge  radius  should  be  as  smell  as 
possible  to  reduce  the  nose  drag. 

To  re^ce  the  required  leading  edge  radius, 
Nonweller**  has  proposed  adding  conducting 
material  aft  of  the  leading  edge  to  transport 
thermal  energy  awey  fron  the  region  of  high  con¬ 
vective  heating  near  the  stagnation  or  attachment 
line,  and  conduct  it  downstream  to  areas  where 
convective  heating  is  lover,  and  excess  energy  can 
be  radiated  away  from  the  body,  honweiler  labeled 
this  theoretical  concept  a;  a  ’conducting  plate’, 
which  is  somewhat  analogous  to  other  passive 
cooling  techniques,  such  as  heat  pipes.  Using 
Nonweiler's  basic  technique,  minimum  leading  edge 
radii  can  be  ascertained,  once  flight  Mach  number, 
freestream  conditions,  leading  edge  sweep, 
material  properties,  and  maximum  allowable  tem¬ 
perature  are  known. 


In  the  present  work,  Nonweiler's  technique 
was  used  to  determine  the  leading  edge  radii  for 
waveriders  designed  for  Mach  numbers  between  6 
and  F*.  The  leading  edge  material  used  for  the 
calculations  was  ATj  graphite,  chosen  because  It 
Is  representative  of  materials  with  high  conduc 
tivity  and  high  melting  point  temperature.  Details 
of  tnis  technique  as  applied  to  the  present  work 
are  given  In  Ref.  7.  It  Is  Interesting  to  note 
that,  for  conditions  associated  with  the  typical 
flight  path  of  a  lifting  hypersonic  vehicle  enter¬ 
ing  the  earth's  atmosphere  at  Mach  25  and  dece¬ 
lerating  to  Mach  6  at  lower  altitude,  the  minimum 
leading  edge  diameters  ranged  from  6  to  28  m*  ~ 
quite  steal  I  in  comparison  to  a  typical  overall 
length  of,  say  60m.  Therefore,  the  present  waveri¬ 
ders  are  essentially  “aerodynamical ly  sharp*  from 
that  point  of  view.  Regardless  of  the  apparently 
small  amount  of  required  leading  edge  bluntness 
(from  the  aerodynamic  heating  point  of  view),  the 
present  waveriders  were  geometrically  altered  to 
accomodate  the  blunt  edge,  and  the  contribution  to 
aerodynamic  forces  on  the  waverider  were  estimated 
assuming  a  modified  Newtonian  pressure  distribution 
on  the  leading  edge. 

0.  Viscous  Flow  Analysis 

A  major  aspect  of  the  present  investigation 
is  that  optimum  waverider  shapes  are  obtained 
wherein  detailed  viscous  effects  are  included 
within  the  optimization  process  Itself.  These 
viscous  effects  are  calculated  by  means  of  two 
Integral  boundary  layer  techniques,  described 
below.  In  all  cases,  the  boundary  layer  flow  Is 
assumed  to  be  locally  two-dimensional,  following 
tha  inviseld  upper  surface  and  lower  surface 
streemlines.  Both  laminar  and  turbulent  flow  are 
considered,  along  with  a  transition  region  based 
pn  empirical  corrections. 

0 . I  Laminar  Analysis 

The  laminar  boundary  layer  calculations 
were  performed  using  Walz  '  integral  method,  as 
described  in  Ref.  24.  The  method  requires  the 
solution  of  a  set  of  coupled  first-order  ordinary 
differential  equations  along  the  boundary  layer 
edge  streamlines.  These  equations  are  the  boun¬ 
dary  layer  moinentum  and  mechanical  energy 
equations,  given  by 


Momentum: 
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Note  that  In  the  above  equations,  primes  denote 
dl fferentlation  with  respect  to  x,  here  repre¬ 
senting  the  boundary  layer  coordinate  In  the 
streamline  direction.  The  variables  In  £qs. 
(9)-(12)  are  defined  by  Walt  to  be 
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Caraful  inspection  of  Eq.  (21)  reveals  that  It  Is 
an  Implicit  equation  for  W*.  since  *  Is  a  complex 
function  of  W*.  Therefore,  In  practice,  a  numeri¬ 
cal  zero-finding  routine  Is  used  on  Eq.  (21)  to 
find  the  value  of  w*  that  yields  the  known  value  of 
W.  w«iz  suggested  an  approximation  that  would 
allow  closed  form  solution  of  Eq.  (21),  however, 
the  present  authors  have  found  that  It  made  boun¬ 
dary  layer  calculations  blow  up  when  used  for  a 
Mach  six  flat  plate  test  case.  Hence,  the 
suggested  approximation  was  discarded  in  favor  of 
the  zero-finding  approach.  For  more  details  on  the 
numerical  solution  of  these  integral  boundary  layer 
equations,  see  Ref.  7. 

0.2  Turbulent  Analysis 

If  and  when  boundary  layer  transition  Is 
predicted,  turbulent  bouhdary  layer  calculations 
are  performed  using  the  Inner  variable  Integral 
method  of  Hhlte  and  Christoph  .  in  practice,  the 
method  requires  solution  of  one  of  two  first-order 
ordinary  differential  equations  along  the  boundary 
layer  edge  streamlines,  depending  upon  the  value  of 
the  parameter  x/*^,  where 
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In  Eq.  (31),  the  parameters  A  and  B  are  defined 
as 
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and  Taw  Is  defined  by  Lq.  (20),  except  that  now 
the  recovery  factor  is  the  turbulent  value. 
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According  to  A* f.  24,  If  A/A  <  0.16,  or  Re*<0, 
then  the  differential  equation* 
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applies,  where 

f*  *  (2.434z  ♦  l.443zZ)exp(-44  z6)  (41) 

g*  ■  1  -  2.3z  ♦  1.762*  (42) 
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and  the  primes  denote,  as  In  the  laminar  case, 
differentiation  with  respect  to  the  streamline 
coordinate,  x.  For  more  details  concerning  the 
numerical  solution  of  tnese  equations,  see  again 
Ref.  7. 


where  x»e  and  xtj  are  the  distances  along  a 
streamline  from  the  leading  edge  to  the  beginning 
and  end  of  transition,  respectively,  and  (Rex)tt 
the  local  Reynolds  number  at  the  beginning  of  tran¬ 
sition  obtained  in  the  present  analysis  from  Eq44  5), 

The  variation  of  local  Skin  friction  coef¬ 
ficient  within  the  transitional  region  (between  *te 
and  xtj)  is  assumed  to  be  a  linear  combination  of 
the  laminar  (cpL)  and  turbulent  ( cy T )  values  that 

would  have  existed  if  the  boundary  layer  were 
completely  laminar  or  turbulent,  respectively.  The 
transitional  friction  coefficient,  Cfra,  is  thus 
related  to  cfL  and  cyT  by: 

S*  •  U-c)  *fL  ♦  e  cfT  (47) 

where  c  is  a  weighting  factor  (a  function  of  x) 
inspired  by  Emmons  (as  discussed  in  Ref. 29).  For 
the  present  Investigation,  the  following  expres¬ 
sion  for  c  is,  as  derived  in  detail  In  Ref.  7: 

<(.,  .  (.8, 

It  is  not  possible  within  the  current  state-of- 
the-art  to  evaluate  the  accuracy  of  these  tran¬ 
sition  correlations.  After  a  study  of  the  existing 
literature,  the  present  authors  feel  that  the  above 
relations  form  a  practical  method  for  simulating 
transition  within  the  goals  of  the  present  study. 
They  provide  a  mechanism  for  assessing  the  effect 
of  transition  on  optimum  waverlder  shapes;  indeed, 
as  discussed  In  the  results,  one  series  of  numerl- 
cal  experiments  Is  conducted  wherein  the  transition 
location  Is  varied  as  a  parameter. 


0.  3  Transition  Analysis 


£.  Aerodynamic  Forces 


The  prediction  of  transition  from  lami¬ 
nar  to  turbulent  flow  at  hypersonic  speeds  Is  a 
state-of-the-art  research  topic.  In  the  present 
analysis,  the  correlation  used  for  predicting  the 
onset  of  transition  is  based  on  two  sets  of  data: 
(1)  data  for  sharp  cones  at  zero  angle-of- 
attack  ;  and  (2)  data  for  win«  with  blunt, 
swept  supersonic  leading  edges7  .  The  correla¬ 
tion  gives  local  transition  Reynolds  number  Rex., 
as  a  function  of  local  edge  Mach  number,  Mj,  as1 
follows  : 

logi0(Rext)  »  6.421  exp  ( l.?09xl0'V  •l"1 )  (44) 


The  lift,  drag,  and  hence  I/O  Is  calculated 
from  a  detailed  integration  of  the  local  surface 
pressures  and  shear  stress  over  the  waverlder 
surface.  Consistent  with  wind  tunnel  practice 
as  well  as  other  literature,  base  drag  Is  not 
included  In  the  present  results,  (for  example, 
all  the  data  shown  in  Fig.  3  do  not  include 
base  drag.)  This  is  done  to  enable  a  rational 
comparison  with  other  data.  Moreover,  at  very 
high  Mach  number,  the  base  drag  becomes  a  small 
quantity  In  comparison  to  forebody  drag.  Details 
on  the  pressure  and  shear  stress  numerical 
integration  can  be  found  In  Ref.  7. 


In  turn,  this  value  of  transition  Reynolds  number 
is  modified  for  wing  leading-edge  sweep,  as 
follows  : 
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where  A  Is  the  sweep  angle,  and  (ReXt)A«0  Is 
obtained  from  Eq.(44).  1 


Once  the  onset  of  transition  has  been  pre¬ 
dicted,  the  extent  of  the  transition  region, 
hence  the  end  of  transition,  Is  predicted  using  a 
relationship  developed  by  Harris  and  Blanchard71 , 
as  follows:  „  , 

(46) 
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F.  Waverlder  Optimization 

Once  a  specific  shape  for  the  forward  leading 
edge  projection  of  a  waverlder  Is  chosen,  (such  as 
shown  in  Fig.  7),  the  techniques  outlined  In  the 
previous  sections  can  be  used  to  generate  the  cor¬ 
responding  waverlder  and  evaluate  its  Hft-to-drag 
ratio  (l/D).  Finding  the  leading  edge  shape  that 
maximizes  l/D,  with  all  other  parameters  held 
fixed,  then  requires  an  optimization  scheme  that 
can  systematically  change  the  projected  leading 
edge  shape  in  search  of  the  one  that  yields  maximum 
l/D.  Unfortunately,  most  existing  optimization 
schemes  require  that  the  function  of  interest  have 
an  analytical  description  —  a  requirement  not  pos¬ 
sible  in  the  present  work.  There  Is  one  scheme, 
however,  a  non-linear  simplex  method  for  function 
minimization  by  Nelder  and  Mead1  ,  that  requires 
nothing  more  than  the  ability  to  numerically 
evaluate  the 
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function.  This  scheme  has  been  used  In  the  p re¬ 
sent  work  to  find  optimum  waveriders. 

In  general ,  the  scheme  of  Ref. 17  minimizes  a 
function  of  n  variables  by  comparing  values  of 
the  function  at  (n*l)  vertices  of  a  "simplex", 
then  replacing  the  vertex  with  the  highest  func¬ 
tion  value  by  another  point  determined  via  the 
logic  of  the  scheme.  As  a  result  of  the 
algorithm  logic,  “the  simplex  adapts  itself  to 
the  local  landscape  [of  the  function  surface], 
elongating  down  long  inclined  planes,  changing 
direction  on  encountering  a  valley  at  an  angle, 
and  contracting  in  the  neighborhood  of  a 
minimum",  according  to  Ref. 17.  In  this  scheme, 
three  operations  —  reflection,  contraction  and 
expansion  —  are  used  to  modify  the  current 
simplex  In  an  attempt  to  replace  the  vertex 
having  the  highest  function  value  with  one  having 
a  lower  value.  Each  of  three  operations  replace 
one  or  more  of  the  (n*l)  points  (Po,Pi .... ,Pn) 
that  define  the  current  simplex  In  n-oimenslonal 
space  with  new  points  that  yield  progressively 
smaller  function  values  (f«,fi»...,fn)  at  the  new 
vertex  points.  A  graphic  illustration  of  how  the 
method  works  Is  shown  In  Fig.  9  for  a  hypotheti¬ 
cal  function,  f,  of  two  variables,  C}  and  C2.  In 
the  figure,  a  triangle  with  vertices  on  the  func¬ 
tion  surface  represents  a  possible  simplex.  In 
the  optimization  process,  the  triangle  (simplex) 
flip-flops  down  the  function  valley,  expanding  if 
possible  to  speed  up  the  process,  then  contract¬ 
ing  when  It  straddles  the  minimum. 

To  use  the  simplex  method  for  optimizing 
waverlder  I/O,  the  shape  of  the  forward  projec¬ 
tion  of  the  leading  edge  must  be  parameterized  In 
some  general  way.  In  the  present  work,  five 
points  In  the  x-y  plane,  lying  Inside  of  the 
shock  domain,  were  chosen  to  represent  the  for¬ 
ward  leading  edge  projection.  A  cubic  spline-fit 
through  the  five  points  Is  then  used  to  generate 
a  continuous  leading  edge.  One  of  the  five 
points,  the  symmetry  plane  point.  Is  constrained 
to  lie  on  the  y-axts,  hence  its  x-  value  Is 
always  zero.  This  leaves  nine  variables,  the 
remaining  x  and  y  values  of  the  leading  edge  pro¬ 
jection  points,  for  the  optimization  routine  to 
manipulate  In  search  of  an  optimum  waverlder,  A 
set  of  leading  edge  coordinates  thus  represent  a 
single  vertex  point, 

Pi  ■  (*2.*3.*4.*5.y1.y2.y3.y4.y5)1  M9) 

of  the  required  simplex,  where  xj  ■  0  as 
explained,  and  the  function  to  be  minimized  Is 
the  negative  of  the  llft-to-drag  ratio 

VV  •  (-1/0),  (50) 

Note  that  the  five  leading  edge  points  are  used 
to  define  only  half  of  the  projected  leading  edge 
shape,  since  the  other  half  Is  constrained  by 
vehicle  symmetry  to  be  the  mirror  Image  of  the 
first  half. 

With  nine  variables  (n-9),  ten  points  (hence 
ten  leading  edge  shapes)  must  be  chosen  to  create 
the  Initial  simplex.  In  the  present  work,  five 
polynomials  of  the  form 

*ze  "  Cl  4  C2*ze  4  C3X*  («> 


yt,  •  C4  4  c5cot(C6  }^li)  (52) 

were  used  to  describe  the  Initial  leading  edge 
shapes;  the  constants  Cj,C2,...C6  being  varied  to 
generate  a  set  of  distinct  shapes.  An  exan^le  of  a 
set  of  Initial  leading  edge  shapes  is  shown  tn  Fig. 
10  —  the  bold  line  representing  the  final  shape 
associated  with  the  optimum  wavender  for  this 
case.  Also  note  that  in  the  present  work,  100 
steps  of  the  optimization  routine  were  executed  for 
all  cases  run,  though  a  convergence  criterion  could 
have  been  Implemented  as  described  in  Ref.  17.  It 
was  found  that  one-hundred  steps  provided  adequate 
convergence  for  engineering  accuracy  ('  10  -lO"** ) 
without  using  excessive  computer  resources  to 
generate  an  optimized  waverider. 

For  more  details  on  the  optimization  scheme, 
see  Ref.  7. 

III.  RESULTS  AND  DISCUSSION 

The  present  results  are  divided  Into  four  sec¬ 
tions,  as  follows:  (1)  a  presentation  of  optimum 
waverlder  shapes  and  aerodynamic  characteristics  at 
Mach  6  and  25,  representing  two  extremes  of  the 
hypersonic  flight  spectrum;  (2)  a  numerical  experi¬ 
ment  to  assess  the  impact  of  boundary  layer  tran¬ 
sition  on  the  optimized  waverlder  shapes;  (3)  an 
assessment  of  the  need  to  account  for  detailed  sur¬ 
face  variations  of  shear  stress  in  contrast  to  the 
use  of  an  average  skin  friction  coefficient  during 
the  optimization  process;  (4)  an  examination  of  the 
question:  If  the  skin  friction  is  deleted  from  the 
present  analysis,  what  type  of  optimized  Invlscld 
waverlder  configuration  Is  produced? 

Due  to  the  specialized  nature  of  my  waverlder 
generation  analysis.  Including  the  present  one.  It 
Is  difficult  to  obtain  a  direct  benchmark  com¬ 
parison  with  existing  data  In  order  to  verify  the 
Integrity  of  the  current  result*.  However,  with 
the  present  analysis.  It  Is  possible  to  calculate 
the  aerodynamic  properties  of  a  half-cone  with  a 
flat  delta  wing  mounted  on  top;  In  this  case  the 
wing  will  have  a  sweep  angle  corresponding  to  the 
shock  angle  of  the  cone,  and  the  body  will  be  at 
zero  degrees  angle  of  attack.  This  specialized 
case  was  calculated  at  Mach  6.8  for  a  half-cone  of 
ec  »  3.67*,  and  the  corresponding  wing  sweep  angle 
of  81*.  The  result  Is  given  as  the  flagged  solid 
square  In  Fig. 3.  This  Is  to  be  partly  compared 
with  the  point  labeled  P2a,  which  was  obtained  from 
Ref.  3,  and  which  corresponds  to  a  similar  flat- 
top  half-cone,  delta  wing  model,  but  at  conditions 
of  maximum  L/0,  hence  at  some  positive  angle  of 
attack.  About  the  only  point  to  be  made  here  Is 
that  the  calculated  L/0  at  zero  angle  of  attack  Is 
lower  than  the  measured  U/0)M.  at  some  angle-of* 
attack  —  a  proper  qualitative^ result.  The 
measured  L/0  at  zero  angl e-of -attack  Is  not  pre¬ 
sented  In  Ref.  3;  however,  through  a  personal 
inquiry  to  Patrick  Johnston  at  NASA  Langley,  the 
present  authors  have  been  told  that  the  measured 
L/0  at  zero-angle-of-attack  was  2.7  —  about  eight 
percent  higher  than  the  value  of  2.5  calculated 
with  the  present  analysis.  This  Is  a  reasonable 
comparison,  and  If  anything,  seems  to  indicate  that 
the  present  aerodynamic  analysis  Is  conservative. 
(Please  note  that  the  comparisons  discussed  above 
are  for  a  given  configuration,  not  an  optimized 
waverlder;  hence  any  degree  of 
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of  validation  here  pertains  to  the  aerodynamic 
portion  of  the  analysis  and  not  to  the  present 
optimization  process  Itself.) 

A.  Representative  Waverlders 

tn  Ref.  7,  a  series  of  optimized  waverlders  Is 
generated.  Including  cases  at  N»  ■  4,6,10,15,20 
and  ?S.  The  conditions  correspond  to  altitude- 
velocity  points  along  a  typical  entry  flight  tra¬ 
jectory  of  a  lifting  hypersonic  vehicle,  such  as 
an  aerospace  plane.  In  the  present  section,  only 
the  results  at  •  6  and  25  are  presented  as 
representative  of  the  two  extremes  of  the  flight 
spectrum.  Ref.  1  should  be  consulted  for  addi¬ 
tional  results. 

Fig.  U  gives  values  of  (1/0),  C^,  and  volu¬ 
metric  efficiency,  n  »  v2/ 3/Sp .  for  waveriders 
optimized  at  different  assumed  wave  angles  for  the 
conical  shock.  To  understand  this  more  fully, 
consider  the  conical  flow  field  associated  with  a 
given  conical  shock  wave,  say  8S  •  11*.  For  this 
value  of  8.  an  optimum  waverider  shape  Is  obtained 
(refer  again  to  the  bold  curve  In  Fig.  10).  The 
resulting  characteristics  of  this  optimized  wave- 
rider  are  then  plotted  on  Fig.  11  for  as  ■  11*. 
This  process  is  repeated  for  other  values  of  8S, 
say  12® ,  13*,  and  14*.  For  each  value  of  8S,  an 
optimized  waverider  Is  obtained,  and  its  charac¬ 
teristics  plotted  In  Fig.  11  as  the  open  symbols. 
(The  solid  symbols  will  be  discussed  later.) 

Hence,  Fig.  il  pertains  to  an  entire  series  of 
optimized  waverlders.  However,  note  that  the  (1/0) 
curve  Itself  has  a  maximum  (In  this  case  for  •«  « 
12*).  This  yields  an  ‘optimum  of  the  optimums*, 
and  defines  the  final  viscous  optimized  waverider 
at  M»  «  6  for  the  flight  conditions  shown  in  Fig. 
11.  The  front  views  of  the  opt  1mm  shapes  at  each 
value  of  ss  are  shown  In  Fig.  12,  and  the  corres¬ 
ponding  perspective  views  are  shown  in  Fig.  13. 
Finally,  a  summary  three-view  of  the  best  optimum 
(the  “optimum  of  the  optimum*)  waverider,  which 
here  corresponds  to  8S  ■  12*.  Is  given  In  Fig.  14. 
Also  in  Figs.  12-14,  the  lines  on  the  upper  and 
lower  surfaces  of  the  waverlders  are  Inviscld 
streamlines.  Note  in  these  figures  that  the  shape 
of  the  optimum  waverider  changes  considerably  with 
es.  Moreover,  examining  (for  example)  Figure  14, 
note  the  rather  complex  curvature  of  the  leading 
edge  In  both  the  planform  and  front  views;  the 
optimization  program  Is  shaping  the  waverider  to 
adjust  both  wave  drag  and  skin  friction  drag  so 
that  the  overall  1/0  Is  a  maximum.  Indeed,  It  was 
observed  In  all  of  the  present  results  that  the 
best  optimum  shape  at  any  given  M  results  in  the 
magnitudes  of  wave  drag  and  skin  friction  drag 
being  approximately  the  same,  never  differing  by 
more  than  a  factor  of  two.  For  conical  shock 
angles  below  the  best  optimum  (for  example  8c  * 

11*  In  Figs.  12  and  13),  skin  friction  drag  is 
greater  than  wave  drag;  in  contrast,  for  conical 
shock  angles  above  the  best  optimum  (for  example 
es  »  13*  and  14*  In  Figures  12  and  13),  skin  fric¬ 
tion  drag  is  less  than  wave  drag.  (Note:  For  a 
hypersonic  flat  plate,  using  Newtonian  theory  and 
an  average  skin  friction  coefficient  ,  it  can 
readily  be  shown  that  at  maxlmumum  I/O,  the  wave 
drag  Is  twice  the  friction  drag.) 

The  results  in  Figs.  11-14  pertain  to  ■  6. 
An  analogous  set  of  results  for  the  other  extreme 
of  the  lifting  hypersonic  flight  spectrum  at  ^  • 


25  is  given  In  Figs.  15-18.  The  aerodynamic  char¬ 
acteristics  of  optimum  waverlders  for  e.  •  7*,8*,9* 
and  10*  ‘are  given  as  the  open  symbols  Fig.  15  (the 
solid  symbols  will  oe  discussed  later.)  The 
respective  front  views  are  shown  In  Fig.  16,  and 
perspective  views  In  FI j .  17.  Finally,  the  best 
optimum  Mach  25  waverider  (which  occurs  at  es  ■ 

9°)  Is  sumxarlzed  in  Fig.  18.  Comparing  the  op¬ 
timum  configuration  at  H,  ■  6  (Fig.  14)  with  the 
optimum  conf Iguratlon  at  Mach  25  (Fig.  18),  note 
that  the  Mach  25  shape  has  more  wing  sweep,  and 
pertains  to  a  conical  flowfleld  with  a  smaller  wave 
angle,  both  of  which  are  Intuitively  expected  at 
higher  Mach  number.  However,  note  from  the  flight 
conditions  listed  In  Figs.  11  and  15  that  the  body 
slenderness  ratio  at  •  6  Is  constrained  to  be 
b/t  *  0.06  (analogous  to  a  supersonic  transport 
such  as  the  Concorde)  but  that  b/t  *  0.09  Is  the 
constraint  chosen  at  1  *  25  (analogous  to  a  hydro¬ 
gen  fueled  hypersonic  aeroplane  such  as  the  ftritlsh 
HOTOL).  The  two  different  slenderness  ratios  are 
chosen  on  the  basis  of  reality  for  two  different 
aircraft  with  two  different  missions  at  either  ex¬ 
treme  of  the  hypersonic  flight  spectrum.  Also  note 
In  Figs.  16-18  the  optimization  program  has  sculp¬ 
tured  a  best  optimized  configuration  with  a  spline 
down  the  center  of  the  upper  surface— an  Interest¬ 
ing  and  curious  result,  due  principally  to  the  com¬ 
peting  effects  of  minimizing  pressure  and  skin 
friction  drag,  while  meeting  the  slenderness  ratio 
constraint. 

Return  to  Fig.  15,  and  note  the  solid  symbols. 
These  pertain  to  the  values  of  C^  and  l/D  obtained 
by  setting  the  ratio  of  specific  heats  y  to  l.l  in 
order  to  assess  possible  effects  of  high  tem¬ 
perature  chemically  reacting  flow.  The  solid  sym¬ 
bols  pertain  to  an  optimized  waverider  at  et  •  9* 
with  y  -  l.l.  This  Is  not  necessarily  the  best 
optimum  at  Mach  25  with  y  «  1.1;  rather,  it  is  Just 
a  point  calculation  to  indicate  that  high  tem¬ 
perature  effects  will  most  likely  have  a  signifi¬ 
cant  impact  on  optimized  waverider  generation,  and 
that  such  effects  are  worthy  of  future  Investiga¬ 
tion.  The  detailed  aspects  of  high  temperature 
effects  are  beyond  the  scope  of  the  present  paper; 
additional  discussion  is  given  in  Ref.  7. 

As  a  final  note  In  this  section,  return  to  Fig. 
3,  and  note  that  the  solid  symbols  pertain  to  the 
present  Investigation.  The  flagged  square  has  been 
discussed  earlier  as  the  data  point  for  a  half-cone 
with  a  delta  wing  at  zero  angle-of-attack;  it  Is 
not  part  of  the  present  waverider  family.  The  un¬ 
flagged  solid  square  at  \  ■  4  pertains  to  a  rela¬ 
tively  large  slenderness  ratio  of  0.087,  used  to 
generate  a  waverider  for  wind  tunnel  testing.  The 
remaining  solid  symbols,  the  circles  and  triangles, 
pertain  to  the  present  discussion.  Recall  that  the 
circles  are  for  b/t  •  0.06  (a  Concorde-like  slen¬ 
derness  ratio  for  a  low  Mach  number  configuration), 
and  that  the  triangles  are  for  b/t  «0.09  (a  M0T0L- 
like  slenderness  ratio  for  a  high  Mach  number 
configuration).  In  the  present  section,  we  have 
discussed  results  obtained  at  Mach  6  and  25;  Fig.  3 
shows  these  plus  others  at  Intermediate  Mach  num¬ 
bers.  All  of  these  cases  are  discussed  In  detail 
in  Ref.  7.  However,  in  regard  to  Fig.3,  emphasis  Is 
now  made  that  the  present  viscous  optimized  waverl¬ 
ders  produce  values  of  (L/D)  which  exceed  the  ”1/0 
barrier*  discussed  in  Section  I,  and  shown  as  the 
solid  curve  in  Fig. 3.  Indeed,  the  present  waverider 

I/O  variation  is  more  closely  given  by 
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<L/n>ma*  -nSr11 

shxn  a$  the  clashed  curve  In  Fig.  3.  Note  that 
the  two  points  given  for  M.  ■  20  and  25  deviate 
away  from  the  dashed  curve.  This  is  a  Reynolds 
number  effect.  Recall  that  all  the  Mach  number- 
altltude  points  for  the  present  waverlders  are 
chosen  to  follow  a  typical  lifting  vehicle  flight 
path  through  the  atmosphere.  The  point  at  Mach 
25  is  at  very  high  altitude  (250,000  ft.),  w1t£  a 
corresponding  low  Reynolds  number  (Re  ■  1.4x10  ); 
the  flow  is  completely  laminar.  At  Mach  20,  the 
Reynolds  number  is  12  times  higher,  but  based  on 
the  transition  criterion  discussed  in  Section  II 
the  flow  Is  still  completely  laminar.  Hence,  the 
laminar  skin- friction  coefficient  at  the  Mach  20 
point  In  Fig.  3  is  much  lower  (cf  *  1//Re)  than 
at  Che  Mach  25  point,  with  an  attendant  larger 
(I/O)  at  Mach  20.  In  contrast,  the  point  at 
Mach  15  Is  transitional,  with  regions  of  both 
laminar  and  turbulent  flow,  and  hence  with 
larger  skin-friction  and  a  lower  (L/0).  In 
any  event,  the  results  given  In  Fig.  3  Indicate 
that  the  present  viscous  optimized  waverlders 
produce  high  values  of  (L/0),  and  therefore  are 
worthy  of  additional  consideration  for  hypersonic 
vehicle  application. 

8.  Sensitivity  to  Transition 

Because  the  major  thrust  of  the  present  work 
Is  the  inclusion  of  detailed  viscous  effects  In 
the  waverlder  optimization,  the  question 
naturally  arises:  How  sensitive  are  the  present 
waverlders  to  uncertainties  In  the  location  of 
transition  from  laminar  to  turbulent  flow?  To 
address  this  question,  a  numerical  experiment  Is 
carried  out  wherein  the  transition  location  was 
varied  over  a  wide  latitude,  ranging  from  all 
laminar  flow  on  one  hand,  to  almost  all  turbulent 
flow  on  the  other  hand,  with  various  cases  Inbe- 
tween.  Specific  results  at  Mach  10  are  given  in 
Fig.  19;  here  values  of  (L/0)  are  given  for  opti¬ 
mized  waverlders  as  a  function  of  assumed  tran¬ 
sition  location.  The  point  corresponding  to  the 
transition  correlation  described  in  Section  II, 
0.3,  is  denoted  by  *x*  In  Fig.  19.  Other  points 
In  Fig.  19  labeled  5x,  lOx  and  15x  correspond  to 
transition  locations  that  are  5,  10  and  15  times 
the  value  predicted  by  the  transition  correla¬ 
tion  ,  All  the  data  given  In  Fla.  19  pertain  to 
optimized  waverlders  for  a.  ■  9*.  which  yields 
the  best  optimum  at  Mach  10  for  the  usual  tran¬ 
sition  correlation.  (Note,  however,  that  9t  ■  9* 
may  not  yield  the  best  optimum  for  other  tran¬ 
sition  locations;  this  effect  Is  not  Investigated 
here.)  The  results  In  Fig.  19  demonstrate  a 
major  Increase  In  (L/D)  In  going  from  almost  all 
turbulent  flow  to  all  laminar  flow.  However,  for 
the  case  where  transition  Is  changed  by  a  factor 
of  five,  only  a  21  change  In  L/0  results.  Cven 
for  the  case  where  transition  Is  changed  by  a 
factor  of  ten,  a  relatively  small  change  In  L/0 
of  111  results.  On  the  other  hand,  the  shapes  of 
the  resulting  optimized  waverlders  are  rairiy 
sensitive  to  the  transition  location,  as 
Illustrated  in  Figs.  20  and  Zl.  The  conclusion 
to  be  made  here  Is  that  waverlder  optimization 
Is  Indeed  relatively  sensitive  to  transition 
location,  and  this  underscores  the  need  for 
reliable  predictions  of  transition  at  hypersonic 
speeds. 


C.  On  the  Use  of  Average  Skin  Friction  Coefficients 

The  present  detailed  viscous  analysis  computes 
the  surface  shear  stress  distributions,  and 
integrates  over  the  surface  to  obtain  the  total 
skin  friction  drag.  This  requires  a  substantial 
amount  of  computer  calculations,  and  leads  to  the 
question:  Can  an  overall  average  skin  friction 
coefficient  be  used  within  the  optimization  pro¬ 
cess  rather  than  dealing  with  the  detailed  shear 
stess  distributions?  To  address  this  question, 
consider  the  best  optimum  Mach  six  case  given  In 
Fig.  11,  which  was  orglnally  calculated  with  the 
detailed  shear  stress  distributions.  From  this 
result,  an  average  skin  friction  drag  coefficient 
was  calculated  for  the  complete  configuration. 

Then  the  optimization  code  was  run  again  for  the 
same  Mach  six  case,  now  using  this  average  skin 
friction  drag  coefficient.  The  results  are  given  in 
Fig.  11  as  the  solid  symbols.  Only  a  small  dif¬ 
ference  exists  between  the  two  cases;  indeed,  the 
resulting  waverider  shapes  are  virtually  the  same, 
as  given  In  Ref.  7.  This  Implies  that  if  an 
accurate  average  skin  ^friction  drag  coeWlcient  can 
be  obtained,  the  resulting  optimized  waverlders 
would  be  reasonably  valid.  However,  the  problem 
with  this  method  1$  that  the  Information  needed  to 
obtain  the  average  skin  friction  drag  coefficient 
Is  not  known  aprlorl.  Moreover,  If  other  indepen¬ 
dent  means  are  used  to  obtain  an  approximate 
average  skin  friction  drag  coefficient  anTThls 
approximate  average  value  Is  used  In  the  optimiza¬ 
tion  process,  the  results  can  be  quite  different 
from  those  obtained  from  the  use  of  detailed  shear 
stress  distributions;  see  Ref.  7  for  more 
discussion  on  this  aspect.  This  situation.  In  com¬ 
bination  with  the  sensitivity  to  transition 
demonstrated  in  the  previous  section,  seems  to  dic¬ 
tate  the  necessity  of  using  the  detailed  shear 
stress  distributions  rather  than  some  approximate 
average  value  of  skin  friction  drag  coefficient  for 
obtaining  the  proper  optimized  waverlders. 

0.  Invlscld  Optimized  Waverlders 

As  a  final  note.  It  Is  Interesting  to  pose  the 
question:  if  the  skin  friction  Is  deleted  from  the 
present  analysis,  what  type  of  optimized  Invlscld 
waverlder  configuration,  with  a  constraint  on  slen- 
derness  ratio.  Is  produced?  To  examine  this 
question,  the  present  computer  code  was  run  without 
skin  friction  as  part  of  the  optimization  process, 
covering  the  range  cf  Mach  number  from  6  to  25.  A 
typical  result  for  the  Invlscld  optimized  configu¬ 
ration  Is  shown  In  Fig. 22.  Here  we  see  essentially 
a  wedge-like  caret  waverlder,  such  as  the  classic 
configuration  generated  by  the  two-dimensional  flow 
behind  a  planar  oblique  shock  wave,  as  discussed  by 
Nonweller  ’  .  This  clearly  Indicates  that  the 
optimized  Invlscld  waverlder  with  slenderness  ratio 
as  the  constraint  Is  Indeed  a  caret  wing.  The 
result  shown  In  Fig.  22  Is  produced  by  the  present 
conical  flow  analysis  as  a  'limiting  case",  wherein 
the  optimum  shape  Is  seeking  the  flattest  portion 
of  the  conical  shock  wave.  To  see  this  more 
clearly,  return  to  Fig.  6.  The  resulting  Invlscld 
waverlders  are  being  generated  by  relatively  flat 
streamsurfaees  at  the  extreme  back  and  bottom  of 
the  generating  conical  flow-field  —  where  the 
shock  radius  of  curvature  Is  the  largest  and  the 
flow  Is  closest  to  being  two-dimensional. 
Consequently,  the  Invlscld  configurations  are  tiny 
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shapes  compared  to  the  scale 
of  the  flowfleld  In  Fig.  6,  and  they  are 
'squeezed*  Into  a  tiny  area  at  the  bottom  of  the 
shock  base.  In  turn,  due  to  the  logic  of  the 
existing  conical  flow  code,  only  a  few  pressure 
and  shear  stress  points  are  calculated  on  the 
surface  of  these  tiny  waveriders,  raising 
questions  about  the  numerical  accuracy  of  the 
calculation  of  their  lift  and  drag.  Therefore, 
no  further  discussion  about  the  invlscid  opti¬ 
mized  waveriders  will  be  given  here,  except  to 
emphasize  again  that  a  two-dimensional  caret  wing 
seems  to  be  the  optimum  Invlscid  shape  that  is 
predicted  by  the  present  conical  flow  analysis. 

IV  CONCLUSIONS 
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In  comparison  to  previous  optimized  waverlder 
analyses,  the  present  work  Is  the  first  to 
Include  detailed  viscous  effects  within  the  opti¬ 
mization  process.  From  this  work,  the  following 
major  conclusions  are  made: 

1.  The  resulting  family  of  viscous  hyper¬ 
sonic  waveriders  yields  predicted  high  values 
of  (L/o)  which  break  the  'L/D  barrier" 
discussed  In  Section  I. 

2.  The  optimization  process  for  the  viscous 
waveriders  results  In  distinctly  different 
shapes  compared  to  previous  work  with 
Invlscld-deslgned  waveriders. 

3.  The  fine  details  of  the  viscous  solution, 
such  as  how  the  shear  stress  Is  distributed 
over  the  surface,  and  the  location  of  tran¬ 
sition,  are  crucial  to  the  details  of  the 
resulting  waverlder  geometry. 
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FIG  1:  Newtonian  results  for  a  flat  plate. 
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FIG.  2:  Performance  comparison  of  various 
generic  hypersonic  configurations 
(  from  Ref.  3). 
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FIG  3:  Maximum  llft-to-drag  ratio  comparison 
for  various  hypersonic  configurations. 
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FIG  4:  Schematic  of  a  simple  caret  Ming 
waverlder. 


FIG.  7:  Front  view  of  a  conical  waverlder 
with  the  generating  cone  and  shock. 


FIG.  5:  Spherical  coordinate  system  for  cone 
flow  calculations,  and  cartesian 
coordinate  system  for  waverlder 
design. 


FIG  6:  Perspective  view  of  a  conical 
waverlder  and  generating  shock. 


FIG  8:  Expansion  cylinder  and  expansion 
domain  bounded  by  the  Mach  cone. 
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FIG.  9:  Example  of  the  optimization  process 
for  a  function  of  two  variables. 
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FIG.  10:  Example  of  Initial  and  optimized 
waverlder  leading  edge  shapes. 


FIG.  11:  Results  for  a  series  of  optimized 
waverlders  at  Mach  6. 
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FIG  12:  Front  views  of  a  series  of  optimized  waverlders  at  Mach  6. 
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FIG.  15:  Results  for  a  series  of  optimized 
waverlders  at  Mach  25. 


FIG.  16:  Front  views  of  a  series  of  optlalzed  waverlders  at  Mach  25. 
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FIG.  19:  L1ft-to-drag  comparison  of  optimized  » 
Hach  10  waverlders  designed  with  various 
boundary  layer  transition  criteria. 
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FIG.  20:  Front  views  of  optimized  Mach  10  waverlders  designed  with  various 
boundary  layer  transition  criteria. 
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Abstract 

An  upwind  difference,  finite  volume  algorithm  was  used  to 
numerically  solve  the  full  three-dimensional  Euler  equations  for 
two  flowfields,  Mach  numbers  2.28  and  5.00,  each  containing  a 
helical-type  vortex  impinging  on  an  oblique  shock  wave.  These 
solutions  showed  several  interesting  and  unexpected  results.  In 
neither  case  did  the  jump  across  the  shock  cause  the  vortices  to 
breakdown.  However,  two  ususual  structures  were  observed. 

The  first  was  common  to  both  flowfields  and  was  seen  as  areas 
of  flow  reversal  around  the  outside  of  the  post -shock  vortex. 
The  second  unexpected  structure  was  observed  only  in  the  Mach 
5.00  flowfteld  and  was  seen  as  the  formation  of  a  trough-like 
structure  around  the  outside  of  the  post-shock  vortex.  Further 
study  is  needed  to  determine  the  cause  of  these  structures.  They 
may  possibly  be  a  result  of  vorticity  generation  due  to  the 
nonuniform  pre-shock  flowfield  combined  with  the  shock  curva¬ 
ture.  To  the  best  of  the  author’s  knowledge,  this  is  the  first 
time  these  two  structures  have  been  observed. 

Introduction 

A  recent  study  at  the  University  of  Maryland  has  been  under¬ 
taken  to  numerically  solve  the  interaction  flow  field  between  a 
3-D,  helical-type  vortex  and  an  oblique  shock  wave.  This  paper 
is  based  on  the  work  presented  in  Reference  1 .  The  study  inves¬ 
tigates  Euler  solutions  to  three-dimensional  shock  wave/ 
vortex  interaction  (SW/VI)  at  two  Mach  numbers,  namely 
Mach  2.28  and  5.00.  Figure  1  shows  a  schematic  of  the  flow- 
field.  This  is  the  first  work  known  to  the  authors  in  which  this 
flowfield  has  been  studied. t 


Computational  domain  =  physical  domain 


Fig.  1  Schematic  of  the  flowfield. 
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tThere  is  one  study  in  which  the  interaction  between  a  vor¬ 
tex  and  a  normal  shock  wave  is  examined  (Ref.  2). 


The  study  is  made  up  of  three  parts.  In  the  first  part,  a  nu¬ 
merical  model  is  chosen  which  is  capable  of  simulating  the 
flowfield.  In  the  second  part,  an  input  vortex  is  developed 
which  can  be  fed  into  the  upstream  boundary  of  the  computa¬ 
tional  domain.  These  tsio* parts  are  then  brought  together  in  the 
third  part  of  the  study  where  the  numerical  algorithm  is  used  to 
solve  the  shock  wave  vortex  interaction  flowfield.  Each  of  these 
parts  arc  discussed  in  detail  below. 

As  mentioned  above,  there  has  been  no  work  studying  this 
particular  interaction  problem,  although  Delery  et  al  in  Refer¬ 
ence  2  have  studied  the  interaction  between  an  normal  shock 
wave  and  a  vortex.  They  experimentally  investigated  the  struc¬ 
ture  of  a  helical-type  vortex  in  a  supersonic  flow  and  studied 
vortex  breakdown^  induced  by  a  normal  shock  wave.  They 
also  developed  an  invtscid  numerical  model  to  simulate  the 
results.  Delery  et  al  chose  to  characterize  the  vortex  by  two 
parameters,  the  freestream  Mach  number  and  the  ratio  of  the 
maximum  tangential  velocity  divided  by  the  freestream  axial  ve¬ 
locity.  Figure  2  shows  the  breakdown  process  as  a  function  of 
these  two  parameters. 
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Fig.  2  Vortex  breakdown  curve.  (From  Ref.  2). 


♦t Vortex  breakdown  is  a  phenomena  wherein  a  vortex  drasti¬ 
cally  changes  its  structure  and  takes  on  one  of  several  new  pat¬ 
terns.  The  restructuring  is  characterized  by  the  information  of  a 
stagnation  point  on  the  vortex  axis  (Ref.  3). 
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The  work  contained- in  Reference  2  formed  a  basis  for  the  in¬ 
put  vortex  used  in  our  study.  This  will  be  discussed  below. 
However,  there  is  an  obvious  and  fundamental  difference  be¬ 
tween  our  work  and  Dclery  el  al's  investigations.  In  iheir 
studies  the  vortex  encounters  a  subsonic  flow  alter  passing 
through  the  normal  shock  and  while  in  our  work  the  Irecstream 
How  is  still  supersonic  after  gotng  through  the  oblique  shock 
wave.  This  may  be  why  no  vortex  breakdown  was  observed  for 
the  two  SW/VI  cases  outlined  below. 

There  are  a  variety  of  applications  where  shock  wave/vortex 
interactions  may  take  place.  Figure  3  presents  three  examples. 
The  first  example,  taken  from  Reference  2,  shows  a  vor¬ 
tex/shock  interaction  on  the  leeward  side  of  a  delta  wing.  Mod¬ 
ern  supersonic  aircraft  and  missiles  require  maximum 
maneuverability  resulting  in  large  angles  of  attack.  This,  in  turn, 
results  in  strong  vortices  being  shed  from  fins  and  wings.  These 
vortices  can  have  a  significant  impact  on  the  vehicle  perfor¬ 
mance  and  even  safety  (Ref.  4). 

The  second  example  shown  in  Figure  3,  again  from  Reference 
2,  is  that  of  a  supersonic  inlet  ingesting  a  vortex.  This  is  obvi¬ 
ously  an  undesirable  situation.  A  knowledge  of  the  conse¬ 
quences  of  this  interaction  is  necessary  to  judge  the  severity  of 
the  situation. 


Shock 


al  Transonic  shock/ 
vortex  interaction 
(From  Ref.  21. 


Breakdown 


bl  Supersonic  Inlet 
(From  Ref.  21. 


The  third  example  presented  in  Figure  3  uses  the  interaction 
process  to  enhance  mixing  in  a  supersonic  combustor.  Whereas 
ilic  First  two  examples  arc  potentially  undesirable  situations,  this 
illustrates  a  case  where  the  interaction  between  a  shock  wave 
and  a  vortex  may  have  a  positive  effect.  The  possibility  ol  using 
SW/VI  in  this  wav  was  suggcsicd  io  rhe  firsi  author  by  Dr. 
David  Van  Wic  of  the  Johns  Hopkins  University  Applied  I'hys- 
ics  1  nboratory. 

Discussion 

The  Numerical  Model 

The  flowfield  is  approximated  by  numerically  solving  the  Eul¬ 
er  equations  rather  than  the  full  Navicr-Stokes  equations.  This 
was  done  for  several  reasons:  I)  solving  the  Euler  equations  can 
be  done  at  a  fraction  of  the  cost  in  code  development  and  com¬ 
puter  storage  and  run  time  over  what  is  required  to  solve  the 
full  Navier-Stokes  equations  and  2)  since  very  little  research  was 
found  on  the  SW/VI  process  a  simple  approach  was  desirable. 
A  number  of  references  were  found  in  which  the  Euler  equa¬ 
tions  were  used  successfully  to  model  vortex  dominated  flow- 
fields  (see  References  2,  5,  6).  These  references  and  our  own 
experience  justified  the  use  of  the  inviscid  equations  for  this 
problem. 

After  determining  the  class  of  equations  to  use  it  was  neces¬ 
sary  to  choose  a  numerical  model.  The  scheme  chosen  is  based 
on  the  algorithm  developed  by  Dr.  Peter  Gnoffo  of  the  NASA 
Langley  Research  Center  (Ref.  7).  Gnoffo  describes  the  model 
as: 

“...  a  robust,  finite  volume,  single-level  storage,  implicit 
upwind  differencing  algorithm  which  has  been  documented 
and  tested  on  several  three-dimensional  blunt  body  flows. 

The  algorithm  can  run  at  unlimited  Courant  numbers  but 
requires  the  inversion  of  only  a  5  x  5  matrix  per  com¬ 
putational  cell.”  (Ref.  7) 

Our  experience  showed  that  the  model  was  indeed  robust,  ac¬ 
curate.  and  capable  of  capiuring  shocks  in  a  crisp  manner. 

The  algorithm  is  briefly  discussed  below  for  completeness 
sake.  A  more  detailed  discription  can  be  found  in  Reference  7. 
The  model  formulation  begins  with  the  integral  form  of  the 
Euler  or  inviscid  equations 
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Fig.  3  Shock  wave/vortex  Interaction  applications. 


p 

> 

0U 

puP  +  pT 

r 

pv 

.7  = 

pvP  +  pj 

,  /t  = 

7 

pw 

pwP  +  pF 

F 

pEP  +  pP 

n 

a 

P 


u.  v,  »• 

E 

P 

P 

T,  J,  F 


volume, 
surface  area, 
density, 

velocity  components  in  the  f,  J,  F  directions, 
total  energy  per  unit  mass, 
velocity  vector, 
pressure, 

the  unit  normal  vectors. 


Unlike  finite  difference  methods  which  approximate  the  tlow 
field  by  a  number  a  jind  points,  a  finite  volume  method  divides 
tile  llowticld  up  inlo  cells.  I  qu.tlioii  tit  is  then  apt'hed  to  a  siii 
pic  cell  williin  the  llowlield  resulting  in  llic  following  appiosi 
mate  equation 
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=  indicates  the  n  +■  I  and  n  time  step 
respectively 

=  indicates  the  ith  cell, 

=  the  volume  of  cell  r, 

=  the  lime  step, 

=  indicates  the  mth  face  of  cell  /, 

=  the  number  of  faces  surrounding  celi  /, 
/,♦  l/|  *  • 

=  indicates  the  cell  face  between  cells  i  and 
/'  +  I, 

=  the  surface  area  between  cells  i  and 
i  +  I. 


The  upwind  differencing  is  introduced  in  the  form  of  a  second 
order  dissipation  term  as  follows: 


converged  quickly  to  a  solution  while  doing  a  very  good  job  of 
conserving  mass,  momentum,  and  energy.  And  as  will  be  seen 
below,  ii  also  coincides  lo  a  solution  in  very  adverse  tlow  mum 
lions,  namely,  the  SW.VI  tlow. 

It  should  he  noted  that  lor  cscry  ease  discussed  below  the  so¬ 
lution  demonstrated  good  global  conservation  of  mass, 
momenutum.  and  energy.  All  eases  where  ran  on  the  University 
of  Maryland  Univac  USX)  computer  using  a  flowficld  containing 
45360  cells  (35  cells  in  the  X-dircctiun  and  36  cells  in  the  Y- 
and  /-directions). 

The  Input  Vortex 

The  nature  of  this  study  required  that  a  vortex  be  provided 
as  an  upstream  boundary  to  the  flowfield.  This  input  sortex 
would  then  sweep  downstream  as  the  solution  was  iterated  in 
time.  Therefore,  it  was  necessary  to  find  and/or  develop  a  suit¬ 
able  input  vortex. 

As  mentioned  above, 'there  has  been  surprisingly  little  work 
done  detailing  supersonic  vortices.  Only  one  reference  was 
found  which  developed  a  useable  input  vortex,  namely.  Refer¬ 
ence  2.  The  input  vortex  discussed  in  Reference  2  was  derived 
by  assuming  uniform  axial  velocity,  zero  radial  velocity,  and 
constant  total  enthalpy.  The  tangential  velocity,  Vn  was 
described  in  terms  of  the  classical  Burgers’  vortex: 
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.  ,(2)  =  the  2nd  order  dissipation  term. 
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=  the  radial  direction 
=  vonex  concentration 
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The  second  order  dissipation  term  is  taken  from  Helen  Yee’s 
paper,  “Numerical  Experiments  with  a  Symmetric  High- 
Resolution  Shock-Capturing  Scheme"  (Ref.  8)  and  is  defined  as 
follows: 

</,.  (2)  =  R„  (d>* i  * 
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With  the  above  assumptions  and  Equation  5,  the  radial  momen¬ 
tum  equation  can  be  integrated  for  static  pressure,  thus,  com¬ 
pletely  describing  the  vortex. 

Two  examples  of  the  tangential  velocity  and  pressure  profiles 
resulting  from  the  above  development  are  illustrated  in  Figure  4 
and  compared  with  the  experimental  results  of  Delery  et  al.  No¬ 
tice  that  increasing  K  increases  the  vortex  strength,  as  expected. 


The  superscript  *  means  to  use  the  most  recent  values  available. 
The  diagonal  matrix  L  is  made  up  of  the  eigen  values  of  matrix 
A  (used  to  linearize  dg  with  respect  to  dq),  matrix  R  is  made 
up  of  the  right  eigen  vectors  of  matrix  A,  and  matrix  R ' 1  is 
the  inverse  of  R.  The  expressions  for  these  matrices  are  devel¬ 
oped  using  Roe’s  Property  U  (Ref.  9).  A  detailed  evaluation  of 
them  can  be  found  in  Reference  7. 

The  model  was  validated  by  comparing  its  results  to  analyti¬ 
cal  solutions  of  2-D  fiowfields  such  as  a  flow-field  containing  an 
oblique  shock,  a  slip  line  intersecting  an  oblique  shock,  and  two 
intersecting  shock  waves.  These  fiowfields,  although  simple,  pro¬ 
vide  a  good  test  of  the  algorithm  because  the  abrupt  changes  in 
the  How  pioperties  across  the  discontinuities  typically  cause 
shock  smearing  and  down  stream  oscillations  in  flowfield  prop¬ 
erties. 

Overall,  the  model  results  compared  very  well  with  analytical 
solutions.  The  numerical  solution  displayed  crisp  shocks  with  a 
minimal  amount  of  smear  and  little  downstream  oscillation  of 
the  flow  properties.  The  algorithm  proved  to  be  robust  and 


To  test  this  model,  a  vortex  was  input  into  a  Mach  2.28  free- 
stream  and  allowed  to  sweep  downstream  without  the  interfer¬ 
ence  of  a  shock  wave.  The  solution  reached  steady  state  in 
about  20  time  steps.  The  pressure  and  density  profiles  from  the 
steady  state  solution  are  illustrated  in  Figures  S  and  6.  Other 
property  profiles  followed  a  similar  trend.  The  solution  shows 
the  vortex  changing  with  downstream  distance  indicating  that 
the  simple  vortex  model  was  not  adequate.  Curves  A,  B.  C, 
and  D  represent  profiles  through  the  vonex  at  increasing  dis¬ 
tance  downstream.  This  progressive  weakening  of  the  vonex  was 
characteristic  of  all  the  vonices  tested. 

There  appeared  to  be  a  slowing  of  the  change  in  the  vonex 
with  downstream  distance.  Therefore,  it  was  thought  that  input¬ 
ting  the  final  vonex  profile  from  the  outflow  boundary  as  the 
upstream  boundary  in  a  new  fiowfidd  might  lead  to  a  con¬ 
verged  vonex  (one  which  did  not  change  with  downstream  dis¬ 
tance).  As  shown  in  Figures  7  and  8  this  was  indeed  the  case. 
The  vonex  dearly  evolves  to  a  converged  solution.  It  was  this 
final  vonex  that  was  then  used  as  the  upstream  boundary  for 
the  Mach  2.28  SW/VI  studies. 
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Fig.  4  Comparison  between  the  experimental  vortex  pro¬ 
files  (Ref.  2)  and  the  vortex  model  profilea  for 
Mach  =  2.28  end  angle  of  attack  =  10  degrees. 
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Fig.  5  Vortex  no.  2  density  cross-section  curves 
through  the  vortex  center. 
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Fig.  6  Vortex  no.  2  pressure  cross-section  curves 
through  the  vortex  center. 
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Fig.  7  Vortex  no.  2 A  density  cross-section  curves 
through  the  vortex  center. 
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Vortex  no.  2A  pressure  cross-section  curves 
through  the  vortex  center. 


A  ■.miliar  approach  was  used  in  developing  a  vortex  lor  the 
Mach  5  IX>  .use  Lnlurtunalcly.  unlike  Hie  Mach  2  2H  case,  the 
viirlcs  lies ^ i  lulls  i.iiivereed  to  a  I m.il  prolile  I  igure  slums 
i lie  density  profiles  through  the  vortex  lor  4  svcles  I  ash  cyji. 
lakes  i he  previous  cycle’s  downstream  vortex  prolile  and  uses  u 
at  an  upstream  prolile.  Ihc  vortex  changes  rapidly  at  lirst  bin 
then,  rather  than  coverging,  continues  to  change  at  a  slow, 
steady  rate.  It  was  assumed  that  it  would  eventually  dissipate 
totally  In  spite  ol  this  anomaly,  it  was  decided  to  use  the  vor¬ 
tex  represented  by  curve  I)  in  I  igure  9b  lor  the  Mach  5  (XI 
N\V  VI  .-tsc. 

The  evolution  of  the  simple  vortices  that  were  input  to  the  li- 
nal  vortices  that  were  output  may  be  caused  by  a  three- 
dimensional  relieving  effect.  The  input  vortex  model  suggested 
by  Dclcry  ct  al  (Ref.  2)  is  a  simple  l-D  solution  to  the  Euler 
equations.  It  is  likely  that  when  the  vortex  is  introduced  into  a 
three-dimensional  field,  it  would  weaken  through  the  3-D  coup¬ 
ling  of  the  governing  equations.  Since  the  vortices  did  eventually 
converge  to  a  constant  profile  (or  nearly  constant  profile),  they 
were  deemed  adequate  to  use  as  input  vortices  for  the  SW/VI 
studies. 

The  Shock  Wave/Vortex  Interaction  Results 

Understanding  the  results  of  the  Mach  2.28  and  5.00  SW/VI 
flow  fields  is  difficult  because,  through  the  interaction  process, 


the  II. m  becomes  lulls  three  dimensional  Inherent  in  this  dilTi- 
uiltv  is  trying  to  vivuali/e  the  3  1)  rcsulis  on  a  2  I)  medium. 

I*»  aid  in  the  proses.,  streamline  and  contour  plots  will  be 
used  along  wnh  cross-sections  curves  of  parameters  cut  through 
i lie  contour  plois.  I  ash  contour  plot  can  be  taken  pcrpcndicul.ii 
io  all  three  axes  and  perpendicular  to  the  frccslrcam  flow  be¬ 
hind  the  oblique  shock  wave.  The  cross-section  curves  can  be 
sin  across  ihc  contour  plot  perpendicular  to  each  of  the  two 
uses  It  quickly  becomes  apparent  that  the  examination  of  a 
three  dimensional  Ho" held  "  a  complex  and  'absirious  process. 

keep  ihc  following  in  mind  when  slewing  the  plots:  Each 
contour  plot  will  usually  have  two  inserts  which  will  schemati¬ 
cally  show  the  location  and  orientation  of  ihc  figure  with 
regards  to  the  rest  of  the  flow-field.  Figure  I  shows  ihc  orienta¬ 
tion  of  ihc  axes.  The  contour  and  streamline  plots  will  be 
shown  in  the  Y-Z  plane  only  and  are  viewed  as  if  the  observer 
were  standing  along  the  negative  X-axis  looking  in  the  positive 
X  direction. 

In  discussing  the  flow  field  solution  it  is  necessary  to  under¬ 
stand  the  swirl  direction  of  the  vortices.  The  vortices  in  this 
study  swirl  in  a  counter-clockwise  direction  when  viewed  as  if 
one  were  standing  on  the  negative  X-axis  looking  in  the  positive 
X  direction.  The  left  and  right  side  of  the  vortex  refer  to  the 
vortex  when  viewed  from  the  negative  X-axis  looking  in  the 
positive  X  direction. 


a)  Cycle  1  —  vortex  3  to  3A  b)  Cycle  2  —  vortex  3A  to  38 


Distance  in  the  Z-direction 
c)  Cycle  3  -  vortex  38  to  3C 


Distance  in  the  Z-direction 
d)  Cycle  4  -  vortex  3C  to  3D 


Fig.  9  Density  cross-section  curvet  through  the  vortex 
center  for  4  different  eyelet. 


Shock  Wave/ Vortex  Interaction  Case  No  I 

The  tirst  cave.  designated  S\V  VI  (  asc  No  I.  introduces  the 
vortex  discuwcd  above  into  a  Math  2.2b  t  loss  field  tonlammg  an 
oblique  shock  wave  rcstiltmg  trom  a  wedge  angle  ol  20  degrees, 
futures  10  and  II  show  the  density  contour  plots  at  various  lo¬ 
cations  along  the  X-axis,  f  igures  12  and  13  show  density  cross- 
section  curves  through  the  post -shock  vortex. 

four  observations  can  be  made  based  on  the  density  plots: 

ll  The  vortex  did  not  break  down  alter  passing  through 
the  shock  wave. 

2)  The  shock  surface  develops  a  convex-concave  shape. 

3)  A  fairly  coherent  area  of  high  density  has  developed  lo 
the  left  of  the  vortex  at  the  point  where  the  vortex  pass¬ 
es  through  the  shock  wave.  This  area  continues  down¬ 
stream  along  with  the  vortex. 

4)  A  large  area  of  low  density  developed  on  the  right  side 
of  the  vortex  after  it  passed  through  the  shefek  wave. 

This  area  almost  encompasses  the  vortex  in  the  last 
plane. 

5)  The  vortex  is  deflected  to  the  left  after  passing  through 
the  shock. 

The  pressure  and  Mach  number  plots  are  presented  in  Reference 
1  and  echo  these  observations.  The  cause  of  these  effects  is  due 
to  the  fact  that  the  velocity  component  on  each  side  of  the  vor¬ 
tex  interseas  the  shock  wave  at  different  angles.  This  results  in 
a  shock  curvature  across  the  vortex  and  the  corresponding  dif¬ 
ference  in  the  properties  across  the  shock  on  the  two  sides  of 
the  vortex. 

Figures  14  and  15  show  the  post-shock  axial  contour  plot  and 
the  tangential  velocity  streamline  plot,  respeaively.  Finally, 
Figures  16  and  17  are  cross-section  curves  corresponding  to 
Figures  14  and  15.  Notice  the  asymmetrical  shape  of  the  axial 
and  tangential  profiles.  This  is  to  be  expected,  again,  because 
the  velocity  vector  intersects  the  oblique  shock  wave  at  different 
angles  on  the  right  and  left  sides  of  the  vortex. 

Although,  the  asymmetry  was  expected,  the  appearance  of 
several  streamlines  running  opposite  the  main  vortex  streamlines 
was  not  expected.  These  streamlines  suggest  the  formation  of 
regions  of  reversed  flow  and  can  be  seen  along  the  outside  of 
the  vortex  shown  in  Figure  15.  The  reversed  flow  io  ihe  right  of 
the  main  vortex  can  also  be  seen  in  the  tangential  velocity  pro¬ 
files  as  a  dip  in  the  profile. 

The  cause  of  this  reversed  flow  is  unknown,  however,  it  may 
be  due  to  vorticity  generation  due  to  the  nonuniform  pre-shock 
flowfield  and  the  shock  curvature.  More  study  is  needed  to  con¬ 
firm  this. 

Shock  Wave/Vortex  Interaction  Case  No.  2 

The  second  case,  designated  SW/V1  Case  No.  2,  introduces 
the  above  discussed  Mach  5.00  vortex  into  a  a  flowfield  con¬ 
taining  an  oblique  shock  wave  resulting  from  a  wedge  angle  of 
32.8  degrees. 

This  flowfield  will  be  harder  to  analyze  than  the  previous 
case  for  two  reasons.  The  first  reason  is  because  of  the  smaller 
distance  between  the  shock  and  the  vortex  at  the  outflow  plane 
due  to  the  decreased  angle  between  the  shock  and  the  post¬ 


shock  freexlrcani  flow.  Thus,  the  part  of  the  flowfielJ  of  in¬ 
terest  is  contained  in  a  smaller  area.  Additionally,  the  changes 
in  ilie  properties  across  the  shock  are  greater  than  for  Cave  No 
I  which  tends  to  wash  nut  the  more  subtle  .hanges. 

Paralleling  the  plots  shown  for  the  case  above,  figures  IK 
and  19  show  the  density  contour  plots  at  sarious  locations  along 
the  X-axis,  figures  20  and  21  show  density  cross-scetion  curves 
through  the  post-shock  vortex.  Figures  22  and  23  show  the 
post-shock  axial  contour  plot  and  the  tangential  velocity  stream¬ 
line  plot,  respectively,  figures  24  and  25  are  cross-section  curves 
coi responding  to  figures  22  and  23. 

It  is  obvious  after  looking  at  the  density  plots  that  there  will 
be  a  lack  of  detail  in  the  graphics.  Case  No.  1,  however,  pro¬ 
vides  some  clues  as  to  what  to  look  for.  There  were  some 
similarities  and  one  unusual  difference  between  the  two  solu¬ 
tions.  The  difference  is  a  trough  encompassing  the  main  vortex 
after  it  passes  through  the  shock  which  is  in  addition  to  the 
high  and  low  areas  seen  in  Case  No.  1.  The  following  summa¬ 
rizes  the  observations: 

1)  There  is  a  trough  that  appears  to  encompass  the  main 
vortex  after  the  vortex  passes  through  the  shock. 

2)  Regions  of  reversed  flow  appear  outside  of  the  main 
vortex. 

3)  The  vortex  has  not  broken  down. 

4)  The  shock  curvature,  the  areas  of  high  and  low  proper¬ 
ties,  and  the  vortex  deflection  noted  in  Case  No.  1  are 
again  present. 

The  presence  of  the  trough  in  this  case  and  not  in  Case  No. 

I  is  most  peculiar.  It  may  be  due  to  the  formation  of  a  slip 
surface  developed  while  the  vortex  is  passing  through  the  shock 
wave  or  it  may  be  a  Mach  number  phenomena  caused  by  the 
increased  shock  strength  in  this  case.  It  may  also  be  related  to 
the  regions  of  reversed  flow.  Again,  more  study  is  needed  to 
fully  understand  this  phenomena. 


Conclusion 

This  concludes  the  discussion  of  the  shock  wave/vonex  inter¬ 
action  studies.  The  two  cases  have  shown  a  variety  of  interest¬ 
ing  results.  A  most  important  conclusion  of  this  study  is  that 
the  vortex  did  not  breakdown  in  either  case.  Although  it  was 
not  noted  above,  the  vortex  strength  did  not  change  appreciably 
in  either  case  after  passing  through  the  shock  wave.  This  may 
be  why  neither  vortex  broke  down.  Related  to  this  may  be  the 
fact  that  the  post-shock  freestream  flow  is  supersonic. 

Both  solutions  showed  the  expected  shock  curvature  and  areas 
of  high  and  low  properties  on  the  outside  of  the  main  vortex. 

In  addition  to  these  effects,  in  each  case,  the  vortex  was  deflea- 
ed  to  the  left  as  it  passed  through  the  shock  wave. 

There  were  two  unexpected  and  fascinating  results  of  the 
study.  The  first  one,  common  to  both  solutions,  was  the  forma¬ 
tion  of  regions  of  reversed  flow  after  the  main  vortex  passed 
through  the  shock  wave.  These  structures  seen  in  the  respective 
streamline  plots  and  the  tangential  velocity  profiles  across  the 
post-shock  vortex  are  a  puzzle.  The  cause  of  the  struaures  is 
unknown.  It  is  hypothesised  that  they  may  be  the  result  of  vor¬ 
ticity  generation  due  to  the  nonuniform  pre-shock  flowfield  and 
the  shock  curvature. 


Fig.  1 1  SW/VI  cat*  no.  1  density  contour  plots  In  the  Y-Z  piano  and  at  X -locations  22  to  34. 
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Fig.  12  SW/V1  case  no.  1  density  cross-section  curve  in 
the  Y-Z  plane,  in  the  Z-direction.  through  the 
vortex  center,  and  at  X-location  =26. 


0  5  10  15  20  25  30  35  40 

Dist.  in  the  Z-direction 

Fig.  13  SW/VI  case  no.  1  density  cross-section  curve  in 
the  Y-Z  plane,  in  the  Z-direction.  through  the 
vortex  center,  end  et  X-location  =  34. 


Plot  orientation 

V  Incl'd  planes  ® 


Freestream  flow 
direction 


Contour  from  .55088  to  .930100  contour  interval  of 
10000-001  pt(3,3l  =  .76821 

Fig.  14  SW/VI  case  no.  1  post-shock  axial  velocity  contour  plot  in  the  110  deg.  plane  and 
beginning  at  Y-location  =  15. 


Plot  orientation 
Y  Incl'd  planes 


.  .  .  /  /VV- 

’  /  //  Ss 


>  >  »  .  *  •  i 

i  i 


K  \\  y 


Freestream  flow  direction 


Fig.  15  SWA/1  case  no.  1  post-shock  tangential  velochy  streamline  plot  in  the  110  deg.  plane 
and  beginning  at  Y-location  » 15. 
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Fig.  16  SW/V1  cast  no.  1  post-shock  axial  velocity 

cross-section  curve  in  the  110  dag.  plana  and 
beginning  at  Y-locati on»  IS. 


Fig.  17  SW/V)  case  no.  1  post-shock  tangential  velocity 
cross-section  curves  in  the  110  deg.  plana, 
through  the  vortex  center,  and  beginning  at  Y- 
locatton  »  15. 
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Fig.  18  SW/V1  case  no.  2  density  contour  plots  in  the  Y-Z  plane  and  at  X-k>cations»6  to  18. 
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Fig.  22  SW/VI  case  no.  2  post-shock  axial  velocity  contour  plot  in  the  123  deg.  plane  and 
beginning  at  Y-location  - 10. 
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Fig.  23  SW/VI  case  no.  2  post-shock  tangential  velocity  streamline  plot  in  the  123  deg.  plane 
and  beginning  at  Y-location*  10. 
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Fig.  24  SW/VI  case  no.  2  post-shock  axial  velocity 
cross-section  curve  in  the  123  deg.  plane, 
through  the  vortex  center,  and  beginning  at  Y- 
location  >10. 


Fig.  25  SW/VI  case  no.  2  post-shock  tangential  velocity 
cross-section  curves  in  the  123  deg.  plane, 
through  the  vortex  center,  and  beginning  at  Y- 
iocation  >10. 


The  other  unexpected  result  was  observed  in  the  second  case, 
but  not  in  the  first.  It  was  the  formation  of  a  trough  surround* 
ing  the  main  vortex  after  it  passed  through  the  shock.  This 
phenomena  could  be  related  to  the  above  mentioned  reversed 
flow  regions  or  it  may  be  caused  by  the  formation  of  a  slip-line 
like  surface. 

In  any  event.  :he  surprising  phenomena  observed  as  a  result 
of  the  numerical  solution  to  these  flowfields,  namely  the  areas 
of  reversed  flow  and  the  trough  structure,  have  not  been 
documented  before  to  the  best  of  the  authors'  knowledge.  How¬ 
ever.  caution  must  be  used  in  drawing  conclusions.  This  has 
been  a  preliminary  study  of  a  very  complex  interaction  flowfteld 
using,  by  necessity,  a  fairly  simple  approach.  The  results  must 
be  weighed  against  this  approach.  It  is  hoped  that  these  initial 
findings  will  gererate  further  interest  and  study  which  will,  in 
turn,  either  support  or  challenge  the  results.  Only  through  more 
analytical,  experimental,  and  computational  study  will  this  flow- 
ticld  be  fully  -roerstood. 
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